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[57] ABSTRACT 

In an ellipsometer, a phase-modulated, polarized light 
beam is applied to a sample, electrical signals are ob- 
tained representing the orthogonal planes of polariza- 
tion of the light after it has interacted with the sample 
and the constants of the sample are calculated from the 
two resulting electrical signals. The phase modulation is 
sufficiently small so that the calibration errors are negli- 
gible. For this purpose, the phase modulator phase mod- 
ulates the light within a range of no more than ten de- 
grees modulations peak to peak. The two electrical 
signals are expanded by Fourier analysis and the coeffi- 
cients thereof utilized to calculate psi and delta. 

10 Claims, 2 Drawing Sheets 
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SMALL MODULATION ELLIPSOMETRY 

RIGHTS IN THE UNITED STATES 

GOVERNMENT 5 

This invention was made with Government support 
under contract NAS8-39327 awarded by NASA. The 
Government has certain rights in this invention. 

This application is a continuation of application Ser. 
No. 07/994,197, filed Dec. 12, 1992, now abandoned. 10 

BACKGROUND OF THE INVENTION 

This invention relates to ellipsometry. 

One class of ellipsometers includes a light source, a 
polarizer, a modulator, an analyzer and at least one 15 
intensity detector. In this prior art type of device, the 
light is modulated and sensed by a photocell. The result- 
ing electric signal is used to calculate psi and delta, 
which can be done by a number of known techniques 
including a Fourier expansion of the resulting signal. 20 

In a prior art type of ellipsometer in this class, called 
a large modulation ellipsometer, the modulator varies 
the intensity of the beam over a defined range suffi- 
ciently large to give an easily measured signal. 

This type of ellipsometry has a disadvantage in that 25 
the system must be calibrated to take into account non- 
linear changes in output signals dependent upon light 
intensity or polarization state. For example, the photo- 
cell output signal in some ranges changes nonlinearly 
with respect to changes in the intensity of light and the 30 
light reflected from diffration gratings, if any are in- 
cluded, change nonlinearly with respect to changes in 
the phase. The changes in ranges of intensity and polar- 
ization state are caused by changes in reflected or trans- 
mitted light as the environment of measurement 35 
changes. 

Because of the need for recalibration with changes in 
the environment, corrections may need to be made such 
as in the software to make the necessary corrections or 
adjustment of voltages or adjustment of amplification 40 
level of amplifiers or sensors or attenuation levels of 
signals from photocells or changes in the angle of the 
diffraction grating if one is included or in the polarizer. 
This alignment is generally obtained by measuring the 
psi and delta of a known environment similar to the one 45 
being tested and calibrating the instrument in that 
known environment. 

This disadvantage of prior art type of ellipsometers is 
especially severe when continuous measurements of 
time varying quantities are made because of the need to 50 
calibrate for different surfaces. 

SUMMARY OF THE INVENTION 

Accordingly, it is an object of the invention to pro- 
vide a novel ellipsometry technique. 55 

It is a still further object of the invention to provide a 
novel ellipsometer capable of measuring time varying 
changes in surfaces in real time. 

It is a still further object of the invention to provide a 
novel ellipsometric technique which does not require 60 
calibration when measuring time varying surfaces over 
a wide range. 

It is a still further object of the invention to provide 
an ellipsometric technique and instrument in which the 
calibration constants do not require frequent adjust- 65 
ment. 

It is a still further object of this invention to provide 
an amplitude modulation ellipsometer which uses suffi- 


2 

ciently small modulation parameters so that the calibra- 
tion errors are insignificant. 

In accordance with the above and further objects of 
the invention, an ellipsometer includes a source of light, 
a polarizer, a modulator and a split analyzer. The sam- 
ple receives polarized light from the source and reflects 
or transmits that light with a change in its polarization 
state. The light is modulated with such a small value as 
to cause the calibration errors related to changes in the 
orientation of the polarized light and light intensity 
insignificant. 

A split analyzer: (1) separates the light reflected from 
or transmitted through the sample into two different 
mutually orthogonal components, one of which repre- 
sents light polarized in alignment with the “p” plane 
and the other of which represents light polarized in 
alignment with the “s” plane; and (2) connects each 
component into a different electrical signal. Although 
the two components have mutually orthogonal polar- 
ization, the polarizations are not identical to the “s” and 
“p” polarizations relative to the sample. The separate 
electrical signals from two light beams are utilized to 
determine psi and delta. 

The amount of modulation that renders the calibra- 
tion constants insignificant is usually in the order of 
three or four degrees, but always under ten degrees and 
preferably under 5 degrees of phase modulation. The 
amount can be determined for specific environments 
from a Jones matrix by deciding on the degree of error 
that is tolerable and setting the angles of the sines and 
cosines below that amount for modulation. For the 
purpose of this specification, the amount of modulation 
is considered as the number of degrees of a complete 
wavelength of light. Amplitude modulation may also be 
used to reduce noise by synchronous detection. 

This ellipsometric technique has several advantages, 
such as: (1) it can measure a changing environment, 
such as for example a surface that is in use such as in a 
spacecraft or the like; and (2) there is less time necessary 
for calibrating the instrument to different surfaces. 

DISCRIPTION OF THE DRAWINGS 

The above noted and other features of the invention 
will be better understood from the following detailed 
description when considered with reference to the ac- 
companying drawings, in which: 

FIG. 1 is a block diagram of one embodiment of small 
modulation ellipsometer in accordance with the inven- 
tion; 

FIG. 2 is a fragmentary block diagrammatic view of 
a modification of the ellipsometer of FIG. 1 that forms 
a second embodiment; and 

FIG. 3 is a block diagram of a third embodiment of 
ellipsometer in accordance with the invention. 

DETAILED DESCRIPTION 

In FIG. 1, there is shown one embodiment of small 
modulation ellipsometer 10 called a polarizer-modula- 
tor-sample-analyzer having a light source 12, a polarizer 
14, a modulation system 16, a dual analyzer 20 and a 
harmonic component detector 22. Light from the 
source 12 passes through the polarizer 14 which pro- 
duces a linear polarization state. The polarized light 
passes through the modulation system 16 which alters 
the polarization state of the light at a high ac frequency. 
The modulation amplitude and thus the magnitude of 
the change in polarization state is kept sufficiently small 
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to cause the calibration constants to be negligible as 
determined by the appropriate Jones matrix analysis of 
the system. 

To obtain signals representing the nature of the sur- 
face of a sample 18, the polarized light is caused to 5 
interact with the sample 18 such as being reflected from 
the sample 18 or transmitted through the sample 18. It is 
shown as being reflected from the sample in FIG. 1. 
The polarized light is altered in state by the delta and psi 
of the sample and the altered polarization state is passed 10 
through an analyzer, such as a Wollaston prism split 
analyzer, to separate the orthogonally polarized compo- 
nents of the beam. The orthogonally polarized compo- 
nents of the beam of light are sensed and used to deter- 
mine delta and psi. 15 

Before converting the beams of light to electrical 
signals, the beam of light incident on the sample 18 or 
the beams of light resulting from it at the Wollaston 
prism split analyzer are phase modulated. The modula- 
tion in the preferred embodiment occurs after the the 20 
beam of light passes through the polarizer 14 and before, 
it engages the sample 18 but the beam may be modu- 
lated at any location along its path up to the time it is 
split in the split analyzer or the beams of light may be 
synchronously modulated after being split and before 25 
being converted to electrical signals. For this purpose, 
the modulation system 16 includes a modulator control- 
ler 24 and a quartz rod 26 or other phase modulator 
which is vibrated to delay and advance the light waves 
as they pass through it to the sample 18. 30 

The modulator controller 24 causes vibration of the 
quartz rod 26 at a relatively high frequency but only 
through a phase angle somewhere between a fraction of 
a degree and ten degrees and preferably under five 
degrees. The exact value is chosen to render the effect 35 
of calibration constants caused by nonlinear characteris- 
tics of components such as the photodiodes to be negli- 
gible. This can be determined from a Jones matrix. 

To obtain two electrical signals, each representing a 
different one of the two orthogonal planes of polariza- 40 
tion of the light after leaving the sample, the dual analy- 
zer 20 includes a split analyzer 30 and two sensors 32. 
The split analyzer 30 separates each of the orthogonally 
polarized components of the beam leaving the sample 
into a separate one of two beams. One such beam split- 45 
ter is a Wollaston prism split analyzer 30. The two 
photodetectors 32 each receive only light in its corre- 
sponding orthogonal polarized plane. Each of the pho- 
todetectors shown at 32 is connected to a different ana- 
log to digital converter within the harmonic component 50 
detector 22. 

The harmonic component detector 22 includes the 
computer 34 and an analog to digital converter system 
36 which receives outputs from two conductors, each 
of the two conductors shown as a single line 33 in 55 
FIGS. 1 and 3 being electrically connected to a differ- 
ent one of the photoelectric devices in the dual photo- 
detectors 32. 

proto.cpp 
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The computer 34 computes the Fourier series using 
standard program software and may be any suitable 
type of computer or microprocessor capable of making 
such a calculation. The A/D board 36 includes two 
analog to digital converters, each connected to a differ- 
ent one of the two conductors shown as the line 33. 

A Fourier expansion of the two signals would result 
in a signal of the type indicated by equation 1 or equa- 
tion 2 where the angular frequency is omega, V a is the 
signal from one of the photosensors, V$is a signal from 
the other photosensor, small t is time and each term 
represents one harmonic of the harmonic expansion in 
terms of voltage. 

The expansion may be formed by any of several 
means although in the preferred embodiment a Fourier 
expansion is obtained using a microprocessor and the 
program provided hereinafter. Commercial Fourier 
transform programs that are suitable are available as 
well. Alternate hardware is available for such a Fourier 
expansion such as by filtering the different frequency 
components in separate filters, separating the sine and 
cosine functions for each frequency by phase detection 
and gating the functions to an output conductor or 
using a harmonic coefficient generator in the type de- 
scribed in U.S. Pat. No. 2,752,092. 


^2Va Z Vb\ - Va 0 Vbi - Va { VbQ t 


EQUATION 3 





EQUATION 11 


r This code is written to acquire data with a DSO .data acquisition board 
manufactured by Metrabyte Co. 

The code will acquire data from the modulation ellipsometer and convert it to Psi 
and Delta. 
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For the executable program to be complete and functional at all this file 
(proto.cpp) should be linked with d50.cpp and show.cpp. 

The code is to be compiled with a Borland C++ compiler or compatible 


written by Hassanayn Machlab June-1992 
7 

#include<dos.h> 

#include<complex.h> 

#inc!ude<stdio.h> 

#include<iostream.h> 

#inciude<string.h> 

#include<stdlib.h> 

#inc!ude<conio.h> 

#include<math.h> 

#include’’has.h" 

#inciude<graphics.h> 
extern int D50lnputlnt(int *lnt); 
extern void GetError(); 
extern FILE *D50; 

FILE *out; 

extern grsphinit(); 

extern void DAS50(void); 

extern void plotfourier(float *,int,int,int 1 int 1 int); 

extern complex *getfourierc(float Tawdata.int); 

extern int D500utput(char *Str); 

int fileread(void); 

int directread(fioat *); 

float *GetPsiDelta(int,char); 

void initD50(void); 

void AcquireAndSave(void); 

void insitu(void); 

int datapts; 

int numwaves=1 ; 

float channel0[2030],channel1[2030],DCoffset0,DCoffset1,m=0; 

int data[2030]; 

extern char filename[20]; 

main() 

{ complex *gfcch0,*gfcch1 ; 
char fileout[20],ch; 
int channels=0, choice, j=0; 
initD50(); 
dc: 


// get background noise before acquiring data 

printffBlock the light beam for dc offset measurement. .\n“); 
printffpush any key to take measurement^!"); 
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getch(); 

directread(channelO); // read data from channel 0 

gfcchO=getfourierc(channelO,datapts); // fourier transform the data 
DCoffset0=abs(gfcch0[0]); 
printf(" DCoffsetO=%f\n",DCoffsetO); 

directread(channei1 ); // read data from channel 1 

gfcchl =getfourierc(channel1 .datapts); // fourier transform the data 
DCoffsetl =abs(gfcch1 [0]); 
printff DCoffsetl =%f\n" I DCoffsetl ); 

printf("\nPush any key to continue acquiring dataln"); 
getch(); 

// start acquiring data 
menu: c!rscr(); 

printf("DATA ACQUISITION FROM MEW); 

printffl. acquire data thru interface or read from file\n"); 

printf(”2.acquire several data and save results to a fils\rf); 

printf(”3. acquire one set of data \n"); 

printf("4.display file contents \n"}; 

printf("5. measure dc offset\n"); 

printf("6.insitu\n"); 

printf(”99. exit\n"); 

scanf("%d" 1 &choice); 

reask: 

if(!(choice == 99 || choice == 4)) 

{' 

printffenter the number of cycles to average at (max=100)\n"); 
scanf("%d", &numwaves); 

} 

if(numwaves>100)goto reask; 

switch(choice) 

{ 

easel: 

channels=fileread(); 

GetPsiDelta(channels,’0'); 

goto menu; 
case 2: 

printf("enter filename to store the data in\n”); 

scanf("%s ,, ,fiieout); 

out=fopen(fileout,"w+"); 

AcquireAndSaveQ;; 
fclose(out); 
goto menu; 



9 


5,416,588 


10 


case 3: 

if (D50Output("Set Channels=0\n")) GetError(); 
directread(channelO); 

if (D50Output("Set Channels=1\n")) GetError(); 
directread(channel1 ); 

// plot data 

graphinit(); 

plotfourier(channel0,datapts,25,290,25,215); 
plotfourier(channel1 .datapts, 31 0,600,25,21 5); 
wind ch00(40,290,260,450,0,0,0,0); 
ch00.plotreal(channel0, datapts, "channel-0"); 
wind chi (31 0,600,260,450,0,0,0,0); 
chi .plotrea!(channe!1 .datapts, "channel-1 "); 
getch(); 
ciosegraphQ; 
cnannels=2; 

GetPsiDelta(channeis, , 0'); 

getch(); 

goto menu; 

case 4: 
clrscr(); 

printf("enter filename to display \n"); 
scanf("%s",fiIeout); ’ 1 

if((out=fopen(fileout,"r")) == NULL) 

{ printfferror opening %s\n",fileout); 
getch(); 
goto menu; 

} 

i-i; 

whiie(fscanf(out,"%c",&ch) != EOF) 

{ 

if(ch -= '\n')j++; 
ifO == 25) 

{j=i; 

getch(); // display 25 lines at a time 

} 

printf("%c" t ch); 

} 

fclose(out); 
getchQ; 
goto menu; 

case 5: 
goto dc; 
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case 6: 

printffenter filename to store the data in\n”); 
scanf("%s",fileout); 

out=fopen(fiieout,"w+"); 

insitu(); 
fclose(out); 
goto menu; 

case 99: 

printf("\nprogram terminated by user !!\n"); 
fclose(D50); 
exit(O); 
default: 
goto menu; 

} 

return(O); 

} 

// Read raw data directly from data acquiaition board without using supplied 
software 

int directread(floai ‘channel) 

{ 

int status.i.bits, channels; 
unsigned int seg.off; 
float fsv; 
char cmd[100]; 
rewind(D50); 

if ( D50Output("CIear\n")) GetError(); // Clear the communication 

buffer 

if ( D500utput("stop\n”)) GetErrorQ; 


datapts=numwaves*20; 

sprintf(cmd,"Set Samples=%d\n",datapts+20); 

if (D50Output(cmd)) GetError(); 

// acquire data 


if (D500utput("Acquire \n")) GetError(); II Start the Trace 
do { 

if (D50Output("REad Status\n")) GetError(); //Wait for Trace to be 

done 

if (D50Inputlnt(&status)) GetError(); 
printf("++++\r"); 

}whiie ((status & 0x18) != 0x18); 
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// tranfer data to array 

rf (D500utput("SET ADD=4 \n")) GetError(); 
sprintf(cmd,'TRAnsfer %u %u 
%d\n ,, ,FP_SEG(data),FP_OFF(data) I datapts); 
if (D50Output(cmd)) GetError(); 


// convert binary numbers to equivalent voltages 

if (D500utput("read range\n")) GetError(); 
if (D50lnputlnt(&status)) GetError(); 
switch(status) 

{ case 0: fsv=2.5; bits=pow(2,11); break; //bipolar 
case 1: fsv=5; bits=pow(2,11); break; 
case 2: fsv=10; bits=pow(2,11); break; 
case 3: fsv=5; bits=pow(2,12); break; //unipolar 
case 4: fsv=10; bits=pow(2,12); break; ' 

} 

fcr(i=0;i<datapts;i++) 

{ 

channel[i]=-fsv*data[i]/bits; // convert to voltage 

} 

channels=1 ; 
return channels; 


} 


// initialize data acquisition board 

void initD50() 

{ 

D50 = fopen( "$DAS50 ,, ,"r+"); // open DAS50 for output and 

input 

if ( D500utput("C!ear\n")) GetErrorQ; // Clear the communication 

buffer 

if ( D50Output("stop\n")) GetErrorQ; 

if (D500utput("Set Channe!s=0\n")) GetError(); 

if (D500utput("Set Samples=48\n”)) GetErrorQ; // should be in steps of 16 

if (D50Output("Set Range= +-10v\n”)) GetErrorQ; 

if (D50Output("Set Trigger Mode=5\n”)) GetError();//6- digital negative edge 
if (D50Output("Set Start=After\n”)) GetErrorQ; 
if (D50Output("Set Rate= Int 1e6\n")) GetErrorQ; 

} 
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/y.«****** j ns jtu*** < ^***************** 

// read data continuosly and store to a file 

void insitu(void) 

{ int channels,parameter,numpoints,j; 
unsigned _delay; 
char ans,titie[20]; 
float *P ARAMS; 
reask: 
clrscrQ; 

printf("Disp!ay quantities or graph fq' or 'g’)?\n"); 

printf("enter 'e' to exit\n"); 

ans=getch(); 

printf("enter total number of points to acquire \n"); 
scanf("%d",&numpoints); 

printf("\nenter delay between points in rnilliseconds\n ,! ); 
scanf("%u"&_delay); 

if(ans=='g') 

{ 

param: 

clrscr(); 

printf("0. graph Delta from channelO\n"); 
printf(”1. graph Psi from channelO\n"); 
printf("2. graph Delta from channel 1\n"); 
printf("3. graph Psi from channe!1\n"); 
scanf("%d",&parameter); 

if(parameter == 0) 

{ 

strcpy(title, "Delta-channel 0"); 
if (D50Output("Set Channeis=0\n")) GetEixorQ; 

} 

if(parameter == 1) 

{ 

strcpy(titIe,"Psi-channel 0"); 
if (D500utput("Set Channels=0\n")) GetError(); 

} 

if(parameter == 2) 

{ 

strcpy(title, "Delta-channel 1"); 
if (D500utput("Set Channe!s=1\n")) GetError(); 

} 

if(parameter == 3) 

{ 

strcpy(tit!e,"Psi-channel 1"); 
if (D50Output("Set Channels=1\n")) GetError(); 

} 

if(parameter <0 |J parameter >3)goto param; 
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graphlnit(); 

wind acquire(40,600, 40,450,0, 0,0,0); 

acquire.plotinsitu(numpoints, title, parameter, _deiay); //0->prints Delta 

getch(); 

closegraph(); 

goto end; 

} 

if(ans=='q‘) 

{ 

j=0; 

while(j<=numpoints) 

{ 

if (D50Output(”Set Channels=0\n")) GetError(); 
directread(channelO); 

if (D500utput("Set Channels=1\n")) GetError(); 
directread(channeil); 

channels=2; 

PARAMS=GetPsiDeita(channe!s,’ '); 
printf("%5.2f %5.2f %5.2f %5.2f %6.4f %6.4f %6.4f\n", 
PARAMS[0],PARAMS[1], 
PARAMS[2],PARAMS[3],PARAMS[5],PARAMS[4], 

PARAMS[6j); 

fprintf(out,"%5.2f %5.2f %5.2f %5.2f %6.4f %6.4f %6.4nn'', 
PARAMS[0],PARAMS[1], 

PARAMS[2], PARAMS[3], PARAMS[5], PARAMS[4], 

PARAMS[6]); 

deiay(_delay); 

j++; 

} 

printf("\n Done!!, push any key to continue \n”); 

getchQ; 

goto end; 

} 

jf(ans=’e')goto end; 


if(ans != 'g' &.& ans != ’q' && ans != 'e') 

goto reask; 

end:; 

} 

r reads data from the DA card through directread() and converts them to Psi and 
Delta. Then saves data in desired file. 

*/ 
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void AcquireAndSave(void) 

{ int channels; 

char ch- O'.chtemp; 
float inread, *P ARAMS; 
clrscr(); 
cntemp=ch; 
whiie(inread != -1) 

{ 

printf("\n enter variable reading, enter -1 to quit\n"); 

scanf("%f',&inread); 

if(inread != -1 ) 

{ 

cirscr(); 

printf("enter channel 'O'.'l'or '2' for 0 & 1\n"); 
ch=getch(); 

if(ch 1= 'O' && ch!= '1 ' && ch!= ’2’)ch=chtemp; 

chtemp=ch; 

if(ch == ’O’) 

{ 

if (D500utput("Set Channels=0\n")) GetErrorQ; 
directread(channelO); 

} 

if(ch == *1*) 

{ 

if (D50Output("Set Channe!s=1\n")) GetErrorQ; 
directread(channell); 

} 

if(ch == '2') II read data from both channels 

{ 

if (D500utput("Set Channels=0\n")) GetError(); 
directread(channelO); 

if (D50Output("Set Channels=1\n")) GetError(); 

directread(channe!1 ); 

channels=2; 

} 

printf("channel= %c\n",ch); 

PARAMS=GetPsiDeita(channe!s,ch); // convert to Psi and Deita 
/* *******”****** store data in a 

fprintf(out," %c %5.2f %5.2f %5.2f %5.2f %5.2f %6.4f %6.4f %6.4f\n ,, 1 
ch, inread, PARAMS[0],PARAMS[1], 

PARAMS[2], P ARAMS[3], P ARAMS[5], PARAMS[4J; 

PARAMS[6]); 

} 

} 
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21 

//*********G e tPsiD e it a ’ 

I* This subdirectory converts the raw data to Psi and Delta. 

for this subroutine ChannelO.channell and datapts have to be declared outside 

main 

7 

float *GetPsiDelta(int channels.char ch) 

{ 

complex *gfcch0,*gfcch1; 

double magch0[20],magch1[20],mcalc1 I mcalc2,mcalc; 
double vaO.val 1 va2,va3,va4,vb0,vb1 ,vb2,vb3,vb4- 
float U, Delta, Tau.T au2, PARAMS[1 0],ca, dummy ; 
int signval ,signva2,signvb1,signvb2; 

intj; 

char *p,mp[10]; 
if(ch !- ’) // for insitu 
{ 

// modulation..m. Used in calculations. 
mp[0]=6; 

printf("\n enter modulation (in waves), default - %6.4f radians\t",m); 
p=cgets(mp); //if return is pushed m keeps a default value 
printfCYO; 

if(mp[1] != 0) 

{ 

m=atof(p); 

m=m*2*3. 141 59265; // convert m from waves to rad 

} 

} 

else 

m=.1 *2*3.141 59265; // convert m from waves to rad 

PARAMS[5]=m; 

I* channel 0 7 

/*** ** * — / 

if(channels==2 J| channels==1 && ch=-0') 

{ 

// get the fourier coefficints and put them in gfcchO. 

gfcchO=getfourierc(channelO,datapts)-; - II fourier transform the data 
for (j=0;j<5;j++) 

{ 

magchO0]=abs(gfcchO03); 
if(ch !=") //for insitu 

printf(”va[%d] = %If phase=%lg re=%lg im=%lg\n",j I magchO0] t 

(atan(imag(gfcchO0])/real(gfcchO0])))*18O/M_PI I real(gfcchO0]) I 

•imag(gfcchO[j])); 

} 
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23 

// find the sign for the first and second harmonics 
signva1=+1; signva2=+1; 

if(reai(gfcch0[7])<0)signva1=-1; // sign of unaveraged val coefficient 
if(reai(gfcch0[8])<0)signva2=-1 ; // sign of unaveraged va2 coefficient 

/‘assign values for clarity. These are the final values of the fourier coefficients 
used in the calculations of Psi and Delta 7 

va0=magch0[0]-DCoffset0; 

val =magch0[1 ]*signva1 ; va2=magch0[2]*signva2; 
va3=magch0[3]*(-1 )*signva1 ; va4=magch0[4]*(-1 )*signva2; 

if(ch I- ') // for insitu 
printf("VaO=%If',vaO); 

r 

results for channel 0 

7 

// Calculation of DELTA 

Delta=atan(0.25*va1 *m/va2); 
if(ch !=' ') II for insitu 

if(signva1==-1 && signva2==+1 ) 
printf("\nDe!ta-a= %f deg\t" I (M_PI+Delta)‘1 80/3. 141 59265); 
else 

printf("\nDeita-a= %f degU", Deltas 80/3. 14159265); 

// Calculation of U and Tau 

U=0.5*cos(Deita)*((va0*m*m/va2)-m*m+4); // channel ’a' calculation 
if(U<=-2)U=U*(-1 ); 
if(ch !=’ ') // for insitu 

printf("U= %At",U); 
if(U>=2) 

{ 

Tau=-0.5*sqrt(U 1, 'U-4)+0.5*U; // channel ’a’ calculation 
if(ch !=’ ’) II for insitu 
printf("psi-a= %f\n\n",atan(Tau)*1 80/M_Pl); 

} ‘ 
else 

{ 

Tau=9999; 

if(ch !=”) //for insitu 
printf("Tau is complexVi 11 ); 

} 

if(signva1 ==-1 && signva2==+1)Delta=M_Pl+Delta; 
if(signva1 ==+1 && signva2==-1)Delta=2*M_Pl+Delta; 
if(signva1==-1 && signva2==-1)Delta=M_PI+Delta; 
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PARAMS[0]=Deita*1 80/M_P I; PARAMS[1]=atan(Tau)*180/M_PI; 

} 


/ ~~~~7 

/* channel 1 7 
/ 7 

if(channels == 2 1[ channels==1 && ch=='1') 

{ 

gfcchl =getfourierc(channei1 .datapts); // fourier transform the data 
for G=0;j < 5;j++) 

{ 

magchl 03=abs(gfcch1 01); 
if(ch !- ') // for insitu 

printf("vb[%d] = %If phase=%lg re=%lg im=%!g\n", j, magchl [j], 

(atan(imag(gfcch1 [jl)/rea!(gfcch1 01)))*1 80/M_Pl,real(gfcch1 01), 

imag(gfcchlOD); 

} 

// find the sign for the first and second harmonics 
signvb1=+1; signvb2=+1; 

jf(reai(gfcch1[7])<0)signvb1=-1; II sign of unaveraged vbl 
if(reai(gfcch1 [8])<0)signvb2=-1 ; I! sign of unaveraged vb2 

II assign values for clarity 
vbO=magch1 [0]-DCoffset1 ; 

vb1=magch1[lj*signvb1; vb2=magch1[2]*signvb2; 
vb3=magch1 [3]*(-1 )*signvb1 ; vb4=magch1 [4]*(-1 )*signvb2; 

if(ch I- ') II for insitu 
printf("VbO=%lf',vbO); 


r 

results from channell 
*/ 


// calculation of Delta 


De!ta=atan(0.25*vb1 *m/vb2); 
if(ch !=' *) II for insitu 
if(signvb1==+1 &S. signvb2==-1 ) 
printfnnDe!ta-b= %f deg\t",(M_PH-Delta)*1 80/3. 141 59265); 
else 

printf( ,, \nDeita-b= %f deg\r.Delta*1 80/3.141 5S265): 
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27 

// calculation of U and Tau 

U=0^5*cos(DeIta)*(-{vb0*m*m/vb2)+nn*m-4); // channel ’b’ calculation 
if(U<=-1)U=U*(-1 ); 
if(ch I- ') // for insitu 

printf("U= %f\t",U); 
rf(U>=1 ) 

{ 

Tau=U-sqrt(U*U-1 ); 
if(ch !=' ’) // for insitu 

printf("psi-b= %f\n\n",atan(Tau)*1 80/M_PI); 

} 

else 

{ 

Tau=9999; 

if(ch !=' ') // for insitu 

printf('Tau is complex'^"); 

} 

if(signvb1==+1 && signvb2==-1 )De!ta=M_P!+De!ta; 
if(signvb1==-1 && signvb2==+1)Delta=2*M_Pl+Delta; 
if(signvb1==-1 && signvb2==-1)Delta=M_Pi+De!ta; 

/*-*+-+* k ***** ** a tr t r * M***k* -* *A *+r* ********* ** i trk * 


PARAMS[2]=Delta*180/M_Pi;PARAMS[3]=atan(Tau)*180/M_PI; 

} 

if(channels ==2) 

{ 

// claculation of modulation 

mcalc=8*va2*vb2/(2*va2*vb2-va0*vb2-va2*vb0); 
mcalcl =8*(va2*vb1/(2*va2*vb1 -vaO’Vbl -val *vbO)); 
if(mcalc<0 j| mcalcl <0) 

{ 

if(ch !=' ’) II for insitu 

{ 

if(mcalc<0) 

printf("\nm calculated from v2's & vO’s is a complex 
number, sqrt(%lf)\n",mcalc); 
if(mcalc1 <0) 

printf("\nm calculated from vI's & v2's & vO's is a complex 
number 1 sqrt(%lf)\n",mcaic); 

} 

} 

else 

{ 
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mcalc=sqrt(mcalc); 
if(ch !=*■•) //forinsitu 
{ 

printffm calculated from v2's & vO's m= %if \n”, mcaic/(2*M_PI)); 
printffm calculated from v1's,v2's & vO's m= %lf PI)); 

PARAMS[4]=mcalc; 

PARAMS[6]=mcalc1 ; 

} 

} 


return P ARAMS; 

} 


r the following subroutine reads data from files ganerated by the software 
supplied with the data acquisition board */ 

int filereadQ 

{ 

char a[50]; 

int numsamplesj.m.e.bitsjj.throwpts.i.channels; 

float fsv.b.dummy.max; 

double rate; 

long pos,mant,expon; 

FILE *in; 

reask: 

DAS50(); 

if((in=fopen(fiIename,"r+"))==NULL) 

{ 

printf("error, couldn’t open %s. Push any key to continued-filename)- 

getchQ; 

goto reask; 

} 

a[0]-\0’; 

while(strcmp(a,"SampleSetSize"))fscanf(in,"%s" 1 &a); 
fscanf(in,”%d",&numsamples); 
whlle(strcmp(a,"RateMantissa"))fscanf(ir , ,,"%s M ; . 
fscanf(in I "%ld’ , ,&mant); 



whiIe(strcmp(a,"RateExponent”))fscanf(in,"%s" 1 &a); 
fscanf(in, , ’%ld”,&expon); 
whiIe(strcmp(a, ,, ChanneJs"))fscanf(in,"%s",&a); 
fscanf(in,"%d” Schanneis); 

rate={(double)mant*pow(10,(double)expon))*channels; 

while(strcmp(a,"Polarity M ))fscanf(in I , '%s" l &a); 

fscanf(in,"%s”&a); 

if(!(strcmp(a,"Bipoiar' , )))bits=pow(2 1 1 1 ); 
if(!(strcmp(a > ,, Unipoiar , ')))bits=pow(2, 1 2); 
whiIe(strcmp(a,"FSMantissa”))fscanf(in,"%s” I &a); 
fscanf(in,"%d",&m); 

while(strcmp(a,"FSExponent"))fscanf(in,"'%s",&a); 

fscanf(in,"%d",&e); 

fsv=m*pow(1 0,e); // calculation of full scale voltage 


clrscr(); 

printf("enter number of waves\n"); 


scanf("%d" I &numwaves); 

datapts=(1/(rate*50e3))*numwaves; 

if(numsamples<datapts) 

{ . 

printf("not enough number of samples\n”); 

fclose(in); 

exit(O); 

} 


while(strcmp(a,"Data ,, ))fscanf(in I "%s" l &a); 

graphinit(); 

for(j=0;j<4;j++)fscanf(in,"%f' I &dummy); 

if(channels==1) 

{ 

forQ'=0;j<datapts;j++) 

{ 

fscanf(in,"%f',&channelO0]); 
channelO0]=-fsv*channelOQ]/bits; // convert to voltage 


plotfourier(channelO,datapts, 25,600, 25,215); 
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wind ch0(40, 600, 260, 450,0, 0,0,0); 
ch0.piotreal(channei0,datapts,”channe!-0"); 
} 


if(channeis==2) 

{ 

forQ=0;j<datapts;j++) 

( 

fscanf(in,"%f,%f\n",&channe!0[j] 1 &channel1 Q]); 
channelOQ]=-fsv*channe!O0]/bits; // convert to voltage 
channell []i=-fsv*channe!1 [jj/bits; // convert to voltage 
} 

plotfourier(channel0,datapts,25, 290, 25,215); 
plotfourierjchannell ,datapts,31 0,600,25,21 5); 
wind ch00(40 ) 290,260, 450, 0,0, 0,0); 
ch00.plotreal(channel0,datapts,"chamel-0"); 
wind chi (31 0,600,260,450,0,0,0,0); 
chi .plotreal(channeI1 .datapts, "channel-1"); 

} 

fciose(in); 
getch(); 
ciosegraphQ; 
return channels; 

} 

In FIG. 2, there is a slight modification of the em- In FIG. 3, there is shown another embodiment 10A of 
bodiment of FIG. 1, which is that a quarter wave com- the ellipsometer substantially the same as the embodi- 
pensator 15 is inserted between the first polarizer 14 and ment 10 shown in FIG. 1 except that it incorporates an 
the modulator 16. This produces circularly polarized additional amplitude modulation stage to reduce noise, 
light at the modulator 16 and nearly circularly polarized 45 In embodiment 10A, a source of amplitude modulation 
light at the sample 18. The detector outputs have the signals 46 modulates the source of light 12 A, which may 
same form as in the embodiments of FIG. 1. The delta be a laser, in intensity so as to provide an intensity mod- 
and psi are determined from equations 8, 9 10 and 11. ulated light source to the polarizer 14 rather than un- 

The amplitudes of the dc and harmonic contributions modulated light from the source of light 12 as in the 

may be extracted using analog or digital filters, and 50 embodiment 10 of FIG. 1. The light signal after passing 
lock-in amplification or digital Fourier analysis. Psi and through the dual analyzer is demodulated before being 
delta can be determined in the case of the embodiment applied to the analog to digital converter system 36 by 
of FIG. 1 from equations 3, 4, 5, 6 and 7. In the case of being applied to a mixer 42 which is also connected to 
the embodiment of FIG. 2, equations 8, 9, 10 and 11 may a function generator or to the source 46 to mix the 
be used to determine psi and delta. In these equations M 55 signals and then being applied to a low pass filter 40 in 
is the modulation amplitude. a demodulator 38. This results in a lower noise signal 

It can be determined from these equations that psi and being applied to the converter, 
delta are independent from the calibration errors. As This amplitude modulation of the light source at a 

can be determined from a conventional Jones matrix, psi frequency f c , called the carrier frequency, which is 

and delta are functions of the sine and cosine of the 60 larger than the frequency of polarization modulation, 
modulation angle. Thus, for modulation angles that are when added to the phase modulation, causes all signals 

very small, the sine is so small that the terms containing to be carried at high frequency, thus reducing noise and, 

the sine are negligible and the terms containing the particularly, eliminating dc offset or drift voltages. The 
cosine are substantially equal to one. Because of this, for original construction and operation of the SME remains 
small angles of modulation, the effect caused by mis- 65 as described above in connection with the embodiment 
alignment does not affect the end result in any signifi- of FIG. 1. 

cant way. Generally, the modulation angle which re- In the operation of the embodiment of FIG. 3, the 

suits in independence from alignment errors is less than source is amplitude modulated at high frequency. Am- 

ten degrees and preferably under five degrees. plitude modulation can be achieved by direct modula- 
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tion of a laser source or by use of an external modulator 
such as an acousto-optic or electro-optic modulator. 
The amplitude modulated light passes through the po- 
larizer 14 and polarization modulator 26, reflects from 
or passes through the sample 18, and then passes 5 
through the split polarizer. 

The orthogonally polarized components of the light 
beam are converted to electrical signals by the two 
detectors 32. These signals are composed of the ampli- 
tude modulated carrier signal at frequency f c with side- 10 
bands at the polarization modulation frequency and 
harmonics as well as undesired dc and ac noise and 
interference. 

The electrical signals are demodulated by standard 
demodulation techniques, such as synchronous (coher- 15 
ent) detection implemented as in FIG. 3. Demodulation 
retains only those signals at f c and its sidebands and 
rejects all signals that are far- from f c . This allows the 
detection of the information carrying signal and the 
rejection of the noise and dc offset signals. The signal at 20 
fc is removed by the low pass filter 40. 

Equations 3 to 1 1 are correct to second order in the 
modulation amplitude M. The accuracy of delta and 
tan. psi may be improved arbitrarily by using higher 
harmonics of the two signals Va and Vb. Alternate 25 
expressions such as shown in equations 3 and 4 for M 2 
permit accurate calculation of delta and tan. psi when 
particular terms of the Fourier expansions (e.g. the nth 
term Va«or Vb n ) have a low signal-to-noise. Preferably 
delta and tan. psi are calculated to the fourth order 30 
although better precision is obtained by higher order 
calculations (use of more terms of the Fourier expan- 
sion). 

In the preferred embodiment, the sources 12 or 12A 
are each a 5 milliwatt diode laser sold as model LAS- 35 
200-670-5 by Laser Max, Inc., 207 Tremont St., Roches- 
ter, N.Y. 14608, although a monochrometer with a 
diffraction grating may be used. The polarizer is a glan- 
laser polarizer, model MGT25S5 and the split analyzer 
is a prism polarizer, model MW2A-10-20, both of which 40 
are sold by Karl Lambrecht, 4204 N Lincoln Ave., 
Chicago, 111. 60618. The compensator used in the em- 
bodiment of FIG. 2 is a precision quarter wave plate 
sold under the catalog no. 02WRQ001-670 by Melles 
Griot, 1770 Kettering Street, Irvine, Calif. 927124-6463. 45 
The modulator is a photoelastic phase modulator sys- 
tem sold under model numbers PEM-90 for the system 
and PEM-90C for the controller, I/FS50 for the head 
and with AR coatings by Hinds Corp., 5250 N. E. Elam 
Young Parkway, Hillsboro, Oreg. 97124-6463. The two 50 
detectors are silicon photodiodes, model SD 100-1 3-13- 
022 sold by Silicon Detector Corp., 1240 Avenida 
Acaso, Camarillo, Calif. 93010 and the analog to digital 
converter board is an analog to digital data acquisition 
board for a personal computer sold under the model 55 
DAS-50 by Keithly/Metrabyte, 440 Miles Standish 
Blvd., Taunton, Mass. 02780. 

As can be understood from the above description, 
this ellipsometric technique has several advantages, 
such as: (1) it can measure a changing environment, 60 
such as for example a surface that is in use such as in a 
space craft or the like; and (2) there is less time neces- 
sary for calibrating the instrument to different surfaces. 
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Although a preferred embodiment of the invention 
has been described with some particularity, many modi- 
fications and variations in the invention are possible 
within the light of the above teachings. Therefore, it is 
to be understood that, within the scope of the appended 
claims, the invention may be practiced other than as 
specifically described. 

What is claimed is: 

X. A method of ellipsometry comprising the steps of: 
polarizing a light beam; 
applying the polarized light to a sample; 
developing signals representing orthogonal planes of 
polarization from the light after it has interacted 
with the sample, Whereby two polarized light sig- 
nals are obtained; and 

calculating constants, psi and delta, of the sample 
from the two resulting polarized light signals, 
wherein a phase modulator phase modulates the 
light within a range of no more than ten degrees, 
whereby calibration errors are sufficiently small to 
be disregarded. 

2. A method in accordance with claim 1 in which the 
two signals are expanded by Fourier analysis and the 
coefficients thereof utilized to calculate psi and delta. 

3. A method in accordance with claim 1 in which the 
light beam is amplitude modulated before being trans- 
mitted to the sample and the two signals are demodu- 
lated to provide noise reducing synchronous demodula- 
tion. 

4. A method in accordance with claim 1 in which the 
beam of light is plane polarized. 

5. A method in accordance with claim 1 in which the 
beam of light is circularly polarized. 

6. Apparatus comprising: 
means for polarizing a light beam; 

means for applying the polarized light to a sample; 
means for developing signals representing orthogonal 
planes of polarization from the light after it has 
interacted with the sample, whereby two polarized 
light signals are obtained; and 
means for calculating constants, psi and delta, of the 
sample from the two resulting polarized light sig- 
nals, wherein a phase modulator phase modulates 
the light within a range of no more than ten de- 
grees, whereby calibration errors are sufficiently 
small to be disregarded. 

7. Apparatus in accordance with claim 6 further in- 
cluding means for expanding the two signals by Fourier 
analysis and using the coefficients thereof to calculate 
psi and delta. 

8. Apparatus in accordance with claim 6 further in- 
cluding means for amplitude modulating the beam of 
light before transmitting it to the sample and demodu- 
lating the two signals to provide noise reducing syn- 
chronous demodulation. 

9. Apparatus in accordance with claim 6 further in- 
cluding means for plane polarizing the beam of light. 

10. Apparatus in accordance with claim 6 further 
including means for circularly polarizing the beam of 
light. 
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